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Abstract 



A maximal coupling of two diffusion processes makes two diffusion particles 
Oh I meet as early as possible. We study the uniqueness of maximal couplings under 

a sort of 'reflection structure' which ensures the existence of such couplings. In 

■4— * ■ 

this framework, the uniqueness in the class of Markovian couplings holds for the 
Brownian motion on a Riemannian manifold whereas it fails in more singular cases. 
We also prove that a Kendall-Cranston coupling is maximal under the reflection 
structure. 
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O ■ 1 Introduction 

The concept of coupling is very useful in various problems in probability. Given probability 
measures //i on (£lx,<^i) and \i 2 on (Q 2 ,^ 2 ), we say \i a coupling of H\ and ii 2 , or 
H G ^(pi, H2), when \i is a probability measure on (fix, #i) x (Q 2 , ^2) so that its marginal 
distributions coincide with /ii and \i 2 respectively. That is, fi(Ax x fl 2 ) = ^1(^1) for 
A x e &x and fi(Qx x A 2 ) = ^2(^2) for A 2 e & 2 . 

We consider couplings of a diffusion process ({Z(t)} t > , {F x } x£M ) on a topological 
space X. A coupling P £ ^(F x ,Fy) determines a stochastic process (Zx,Z 2 ) on X x 
X so that each individual component moves as the diffusion process starting at x and 
y respectively. A characteristic of couplings on which we concentrate our attention is 
the coupling time T(Zx,Z 2 ), the time when Z\ and Z 2 coalesce (denned in (12.11) ). In 
many applications, we would like to make the coupling probability P[T > t] small by 
constructing a suitable coupling P. In these ways, one can obtain various estimates for 
heat kernel, harmonic functions(or harmonic maps), eigenvalues etc. by means of the 
geometry of X. These results indicate that the existence of a good coupling reflects the 
nature of 2 or A. 

Our interest in this paper is the problem of the uniqueness. More precisely, we would 
like to know what properties of Z or X are related to the uniqueness of couplings which 
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minimize the coupling probability. At this moment, however, the existence of such a good 
coupling is not obvious at all in general. Thus we confine ourselves in a special situation 
where the existence is ensured. 

In the preceding work by E. P. Hsu and K.-Th. Sturm [Bj, they discussed the uniqueness 
of maximal coupling when X = M. d and Z is the Brownian motion on it. Motivated by 
the coupling inequality, they defined a maximal coupling as it minimizes the coupling 
probability. In their framework, there is a natural maximal coupling Fm £ ^(P^nP^) 
called "mirror coupling" defined by using the reflection with respect to the hyperplane 
which maps X\ to x-i- They showed that the mirror coupling is the unique maximal 
coupling in the class of Markovian couplings ^(P^, P^) (see Definition 12. 4p . They also 
showed by examples that the uniqueness no longer holds when we are allowed to take 
non-Markovian couplings. Their argument to derive the uniqueness uses the explicit form 
of the transition density of the Brownian motion. In this sense, their argument depends 
on the nature of the Euclidean Brownian motion. 

In order to investigate how such a uniqueness depends on the nature of Z or X, we 
discuss the same uniqueness problem in a similar, but more general, situation. That is, we 
assume a sort of 'reflection structure' like a reflection in Euclidean spaces for given initial 
points Xi,X2 € X. Then we can naturally define a mirror coupling Fm £ ^(Pm, F X2 ) as 
a maximal Markovian coupling. In this situation, we consider the uniqueness of maximal 
couplings in c ^'o(F Xl ,F X2 ). As a result, the Brownian motion on a Riemannian manifold 
enjoys the uniqueness. But, as we will see, the uniqueness no longer holds if we consider 
more singular cases. These observations show that the uniqueness is related to the nature 
of Z or X even when the mirror coupling exists. 

The organization of this paper is as follows. In the next section, we introduce our 
framework including the notion of 'reflection structure', maximal coupling and Markovian 
coupling. Our main theorem gives a sufficient condition to the uniqueness of maximal 
couplings in %)(F X1 , F X2 ) (Theorem I2.6p . We will prove Theorem I2.6l in section [3] following 
the idea of [BJ. Section H] is devoted to some examples. On one hand, we will show that 
the uniqueness holds under the assumption on the short time asymptotic behavior of Z 
and the geometry of X (Theorem 14. ip . A typical example satisfying these conditions is the 
Brownian motion on a complete Riemannian manifold (Corollary 14.31) . This framework 
includes the Euclidean Brownian motion as discussed in [6j. There we exhibit complete 
Riemannian manifolds which have the reflection structure with respect to specified initial 
points. On the other hand, we also show two easy examples where the uniqueness of 
maximal Markovian coupling fails (see Example 14.101 and Example 14.111) . At the end of 
this section, we consider the case for the Brownian motion on 2-dimensional Sierpinski 
gasket. We show that the uniqueness holds while this case is not included in the framework 
of Theorem 14. 1[ In section we show that the Kendall- Cranston coupling coincides with 
our mirror coupling under the existence of the reflection structure. The Kendall-Cranston 
coupling is originally introduced by Kendall [8J and Cranston [3j for the Brownian motion 
on an arbitrary complete Riemannian manifold. Their coupling is useful to estimate 
analytic quantities by means of the geometric quantity such as Ricci curvature. But, 
in general, there is no reason why the Kendall-Cranston coupling should be maximal. 
The construction of their coupling is based on a sort of reflection of infinitesimal motion 
by means of the Riemannian geometry. Thus, our result is quite natural. It should 
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be remarked that our result implies that the Kendall-Cranston coupling is the unique 
maximal coupling if there is a reflection structure. 



2 Coupling of diffusions and its properties 

Throughout this paper, we assume X to be an arcwise-connected Hausdorff topological 
space with the second countability axiom. For a coupled diffusion process (Zi(t), Z 2 (t)), 
the coupling time T(Z\, Z 2 ) is defined by 

T(Z U Z 2 ) := mf{t > ; Z x {s) = Z 2 {s) for all s > t}. (2.1) 

We set 

<p t (x, y) : = \ \\F X o Z{t)- X - F y o Z{t)-\ w . (2.2) 



2 

Here ||-|| var stands for the total variation norm. By using this function, the coupling 
inequality is written as follows: for every x, y G X and P G ^(P x , F y ), 

F[T{Z 1 ,Z 2 )>t]>i Pt {x,y). (2.3) 

For the proof of (I2.3p . it suffices to remark that 

F[T(Z 1 ,Z 2 )>t]>F[Z 1 (t)^Z 2 (t)} 

>¥[Z 1 (t)eA,Z 2 (t)^A] 
>¥[Z 1 {t)eA]-¥[Z 2 (t)eA] 

holds for arbitrary A G SS[X\ 

Definition 2.1 (cf. P [6]) For t > 0, we say P G tf(F x ,F y ) maximal at t when the 
equality holds in (I2.3p . We say P G ^(P^jP^) maximal when the equality holds in (12. 3p 
for each t > 0. 

Let us fix Xi,x 2 G X. The reflection structure with respect to x\ and x 2 stated in section 
[T] means the following two properties assigned on X and Z: 

(Al) There is a continuous map R : X — > X with R o R = id so that F Xl o R~ l = F X2 , 

(A2) The set of fixed points H := {x G X ; i?(x) = x} separates X into two disjoint 
open sets X x and X 2 (i.e., X \ H = X 1 U X 2 ) with i?(X x ) = X 2 . 

As an easy but significant consequence of (A2), every continuous path in X joining x G X\ 
and y G X 2 must intersect H. In general, it highly depends on the choice of Xi,x 2 G X 
whether (Al) and (A2) hold or not (see Example 14.81) . But, we can easily verify that (Al) 
and (A2) are satisfied for the Euclidean Brownian motion for any x±,x 2 G X. In that 
case, R is an reflection with respect to a hyperplane H. Under (Al) and (A2), we can 
construct a mirror coupling of F Xl and F X2 . Let r := inf{t > ; Z\{t) G H} be a hitting 
time to H. We define the mirror coupling Fm as the law of (Zi, Z 2 ) where Z\ is a copy 
of (Z,F X1 ) and 

::;>:: 



By definition, F M G tf(F Xl , F X2 ) and r = T(Z U Z 2 ) under P 
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Proposition 2.2 Fm is maximal. 

For the proof, we use the following lemma. 

Lemma 2.3 Suppose F x o R^ 1 = P Ra , for x G X\. Then, for each t > 0, 

<^(x, ifcr) = P x [Z(t) G X x ] - P te [Z(t) G X x \. 

Proof. By Q , 

^(x,i?x) = sup (F x [Z(t)eA]-F Rx [Z(t)eA}). (2.5) 

Ae.«(x) 

Note that 

F x [Z(t) £ A] — F Rx [Z(t) G A] = F x [Z(t) E A, r < t] - F Rx \Z(t) G A, r < t] 

+ F x [Z(t) eA,r>t]- F Rx [Z(t) eA,r> t}. 

First we show 

F x [Z(t) EA,r<t] = F Rx [Z(t) eA,r<t] (2.6) 
for each A G 3B{X). By the strong Markov property, 

F x [Z(t) eA,r<t]= E x [F z(T) [Z(t - s) G A]| S=T ; r < t] . 

By assumption, the law of {Z(t), t) under P^, equals that under P^.. Thus we have 

E,. [F Z (T)[Z(t -s)e A]\ S=T ; r < t] = E Rx [F z{r) [Z(t - s) G A]| S=T ; r < t] 

= F Rx [Z(t) eA,r<t}. 

Next, by (A2), we have 

F x [Z(t) eA,r>t]- F Rx [Z{t) eA,r>t] 

= F x [Z(t) eX 1 nA,T>t]- F Rx [Z(t) eX 2 nA,r> t}. 

These observations imply that the supremum in (12.51) is attained when A — X\. □ 
Proof of Proposition [Ml By (A2), 

F M [T(Z 1 ,Z 2 ) >t]= F Xl [r >t]= F Xl [Z(t) eX 1 ,r>t)- F X2 [Z(t) G X u r > t). 
By (Al), we can apply (12.61) for x = X\. Thus we obtain 

F Xl [Z(t) G Xi, t > t] — F X2 [Z(t) G Xi, t > t] — F xl [Z(t) G X x ) - F X2 [Z(t) G X 1 }. 
Hence Lemma [2.31 yields the conclusion. □ 
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Definition 2.4 Let Z* = (Zi, Z 2 ) be a coupling of diffusion process Z starting from 
(xi,x 2 ) under P. We define a canonical filtration {^*}t>o by := a{Z*(u) ; <u < 
s}. We say that P is Markovian or P G ^ (P X1 ,P X2 ) if, for each s > 0, the shifted process 
{Z*(t + s)}j>o under P conditioned on is still a coupling of the diffusion process 
starting from Z*(s) = (Zi(s), Z 2 (s)). By using the shift operators {0 S } S>0 defined by 
(9 s (Z*))(t) = Z*(s + t), P G %(F Xl ,F X2 ) means P[ • \& a \ oO; 1 G tf(P Zl ( a ), P^w) for each 
s > 0. 

Obviously, the mirror coupling P^ is Markovian. As noted in [5], the condition that Z* 
is a Markovian coupling does not imply that Z* is a Markov process in general. 
To state our main theorem, we introduce a subclass of ^{F X1 , P^). 

Definition 2.5 We say P G tf(F xl , F X2 ) when, for each t > 0, there is E® G ^(X x X) 
with Clixl 2 and P(Z*(t) G (S^) c n X x x X 2 ) = so that each (x,y) G 5 (t) 
satisfies the following: if there is a decreasing sequence {s n }neN of positive numbers with 
lim n ^oo s n = so that 

F x [Z{s n ) eA]> F y [Z{s n ) G A), (2.7) 
F x [Z(s n ) G A'] < P„[Z(s n ) G A'] (2.8) 

hold for all A C XiU H, A' C X 2 U H and all n G N, then x = Ry. 

We can easily verify P M G if (P« x , P^). 

Theorem 2.6 Assume (Al) and (A2)forx 1} x 2 G X. Let P G ^(P^, P^) H^P^, P* 2 )- 
If there is to > so that F is maximal at every t G (0, to); t/ien t/je law of Z*(tAto) under F 
is identical to that under Fm- In particular, ifF G ^(P^, P^) H c &o(F Xl , F X2 ) is maximal, 
then P = Pa/- ^4s a result, if 

^iF Xl ,F X2 ) D < ^,(Px 1 ,Px 2 ), (2.9) 
t/ien Pm t/ie unique maximal coupling in %)(¥ Xl , F X2 ). 

Remark 2.7 (i) The conditions (12. 7p and (12. 8p are equivalent to the fact that a Hahn 
decomposition of F x o Z{s Tl )^ 1 — F y o Z(s n ) _1 is given by Xi or Xi U for each n£N. 
(ii) We can directly show that the Brownian motion on a Euclidean space satisfy (12.91 ). 
Indeed, for every x, y G P d , X = {2 G FL d ; |z — x| < \z — y|} gives a Hahn decomposition 
of Pa; o Z(s)^ 1 — F y o Z{s)~ 1 for each s > 0. This is because the transition density depends 
only on the distance for fixed t > 0. (iii) In the case of Euclidean Brownian motion, more 
strong assertion holds: the maximality of P only at t > implies that the law of Z*(- A t) 
under P is identical to that under F M (see [5]). But their proof requires some properties 
derived from the explicit form of the transition density of the Euclidean Brownian motion. 
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3 Proof of Theorem 12.6 



To begin with, we remark that (Al) produces the following auxiliary lemma. 

Lemma 3.1 For each t > 0, there is X® G 38{X) with F xl [Z(t) G X®] = 1 so that 
P 2 o R- 1 = F Rz for z G X® . 

Proof. Take A f G 38{X) for % — 0, . . . , n and < s 1 < s 2 < ■ ■ ■ < s n . Then the Markov 
property implies 

P X1 [Z(t) G Aq, Z{t + Si) G Ai, . . . , Z{t + s n ) G A n \ 

= E X1 [l Ao (Z(t)) ■ P m [Z( Sl ) G At, . . . , Z{s n ) G A n )} . 

By using (Al) twice, 

F X1 [Z(t) G A ,Z(t + si) G Ai, . . . , Z(t + s n ) G A n ] 

= F X2 [Z(t) G R^Ao, Z(t + Sl ) g iT%, ...,Z(t + s n )e R~ l A n ] 
= E X2 [l R -i Ao {Z{t)) ■ F m [Z{si) G R- l A u Z{s n ) G R~ l A n ]] 
= E X1 [U (Z(t)) ■ F RZ(t) [Z(si) G R~ l A l} . . . , Z{s n ) G iTMJ] . 

Since A is arbitrary, there is X su _ >Sn . Au _ >An G 38{X) with P^ G A Slv .. jSn . Alv .. iA J = 1 
so that 

P, [Z{s x ) G Ai, . . . , Z(s n ) G A n ] = P fe g i?- 1 ^!, . . . , Z(s n ) G R~ l A n ] 

for a; G X Sl) ... tSn;Al) ...,A n - Since A enjoys the second countability axiom, there is a countable 
family of open sets U in X so that cr(U) = 38 (X). Thus 

x v = n n n x sl ,..., sn , Al ,..., An 

neN Si£Q Ai&t 
l<i<n l<i<n 

is what we desired. □ 

Remark 3.2 In this paper, we used the second countability axiom of X only for the 
proof of Lemma [331 Thus, if F x o R^ 1 = P^. holds for all x G A, then X need not satisfy 
it. 

We write fi\ = F Xl o and /4 = P^ ° Zit)~ x for simplicity. Let us define /ig by 

/4(A) = /U *(AnA 1 ) + ^ 1 (AnA 1 c ) (3.1) 

for each A G 38 {X s ). By Lemma [2.31 we have /Xq < /x' and /Jq < /tx|. 

Definition 3.3 For t > 0, the mirror coupling /x^ G ^(/x* , /x|) is the probability measure 
on X x X defined by 

^ M {dxdy) = 5 x {dy)/J {dx) + 5 Rx {dy){fi\ - nl)(dx). (3.2) 
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We can easily verify G - 
Lemma 3.4 Let s,t > 0. T/ien /or x,y E X, 

inf( / Vs {z ll z 2 )v{dz l dz 2 )] ue^iF.oZity^FyoZit)' 1 )) >(p s+t (x,y). (3.3) 
Uxxx J 

In particular, the equality holds when (x,y) = (xi,x 2 ). In this case, the infimum is 
attained at f/ M . 

Proof. Let u tjE (z) : = F z [Z(t) G E] for £ G =^(X). Let a* G ^(P x o Z(t)~ x , F y o Zity 1 ). 
Then 

u s+t , E {x) - u s+t ,E(y) = E x [u s>E (Z(t))] - E y [u StE (Z(t))] 

{u s ,e{zi) - u StE (z 2 )}du t (dz 1 dz 2 ) 



XxX 

tl 



< / ip s (z l7 z 2 )d/j (dzidz 2 ). 



XxX 



By taking the supremum on E G <^(X) in the left hand side of the above inequality, we 
obtain ( 13.31) . We now turn to the latter assertion. We set x = x\ and y = x 2 . By (13. 21) . 
we have 



(p s (z 1 ,z 2 )dn M (dz 1 dz 2 ) = / ip s {z,Rz)^ l {dz) - I ip s (z, Rz)/j Q (dz). (3.4) 

XxX Jl ./X 

Set Ut(z) = ut^iz). Let X® be as in Lemma I37T1 By Lemma l2~3l we obtain 

(p s (z, Rz) = u s {z) - u s (Rz) (3.5) 
for z G X\ fl (XW u itX"W). Thus 

(p s (x, Rx)n t (dx) = / (p s (x, Rx)n t 2 (dx) + / (f s (x, Rx)n t 1 (dx). (3.6) 
x ixi ix 2 

Substituting (13.61) to (13.41) . we obtain 

XxX JXi JXi 

(f s (x, Rx)n\(dx) — / (f s (x, Rx)n t 2 (dx) 

{u s (x) — u s (Rx)} [i\{dx) — I {u s (x) — u s (Rx)} fJ 2 (dx) 
X Jx 1 

u s (x)n t 1 (dx) — / u s (x)f/ 2 (dx) 
x Jx 

= u s+t (xi) - u s+t (x 2 ) 
= <p s+t (xi,x 2 ). 

Here the third equality follows from (13.51) . □ 
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In the following, we show a kind of converse assertion. 

Proposition 3.5 Let P G e ^(F Xl ,¥ X2 ) and t > 0. Suppose that there is a sequence 
{■Sn}neN so that 

E[<p 9n (Z 1 (t),Z 2 (t))] = <p 9n+t (xi,x 2 ) (3.7) 
holds for all n£l Then P o (^(t), Z 2 (*)) _1 = Hm- 

Let Z) := {(x, x) G X x X ; x G X} and i : I-»I)a canonical injection. For the proof 
of Proposition 13.51 we show the following lemma. 

Lemma 3.6 Let P G ^(P^, P X2 ) and t > 0. Suppose that there is a sequence {snjneN so 
that 

H ( Ps n {Z 1 (t),Z 2 {t))} = <p Sn+t (x h x 2 ) 
holds for all n£l Then P o (Zi(t), Z 2 {t))~ x \ D = ^ o t" 1 . 

Proof. Set /i := P o Z 2 (t)) -1 . For simplicity, we write /i' =: /ij for i — 0, 1, 2. By 

a usual argument, /x is expressed in the following forms: 

/i(dxdy) = ki(x,dy)/j l i(dx) = k 2 (y,dx)/i 2 (dy). 

We define a coupling v G ^(/ii, /i 2 ) by 

1 1 /" 

v{dxdy) = -5 x (dy)u. (dx) + - I k 2 (z,dx)ki(z,dy)u. (dz) 
^ ^ Jx 

- -h(x, dy)/j, (dx) - -k 2 (y, dx)u, (dy). + fi(dxdy) 
By (13.31) and ( \3.7L for s G {snjngN, we have 

0< / p s {x,y)v(dxdy) - <p s {x,y)u,(dxdy) 

JXxX J XxX 

= - I <f s (x,y)k 2 (z,dx)k 1 (z,dy)iJ l0 (dz) 



JXxXxX 

~ J <f s (x,y)k 1 (x,dy)fx {dx) - ~J p s (x,y)k 2 (y,dx)u. (dy) 

{ip s (x, y) - cp s (z, y) - (p s (z, x)} k 2 (z, dx)k 1 (z, dy)/j, (dz). (3. 

XxXxX 



By the triangular inequality for ||-|| var , we have ip s (x,y) < ip s (x,z) + <p s (z,y). Thus 
the left hand side of H 3 . 8 1) must be 0. Moreover, there is Q s C X x X x X with 
Jnc k 2 (z, dx)ki(z, dy)fi (dz) = so that 

cp s (x, y) = (f s (x, z) + <p s (z, y) 

holds for each (x,y,z) G f2 s . Note that this equality is equivalent to the existence of a 
Borel subset Ei x ' y ' z ^ C X which satisfies 

V x [Z(t) eA}< ¥ z [Z(t) eA}< ¥ y [Z(t) G A] 
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for each Borel set A C Ei x ' v ' z ^ and 

F y [Z(t) G A'} < F z [Z(t) G A'} < ¥ x [Z(t) G A'}. 

for each Borel set A' C (Es X ' v ' z ^) c . This fact follows from a simple calculation of the total 
variation norm by using Hahn decompositions. Let Q := f] n&N fl Sn . We set 

Ai = {(x, y, z) G X x X x X ; x = z}, 
A 2 = {(x, y, z) G X x X x X ; y = z}. 

Then we claim 

ficAUAs. (3.9) 

Let (zi,Z2,zs) G fi. Suppose (zi,Z2,Zs) ^ U A 2 . Take open neighborhoods Vi of z« 
{i = 1,2,3) with Vi fl V3 = for i = 1,2. We choose n G N sufficiently large so that 
P*i[£(s„) G V<] > 3/4 for i = 1,2,3. But, for E Sn = E { s z n 1,Z2 ' Zs) , we have 

\ < V Z3 [Z{sn) e F 3 ] = P Z3 [Z(s n ) g V3 n £ S J + P 23 [Z(s n ) g V 3 n £ S C J 

< P 22 [z( Sn ) g r 3 n £J + F zl [Z(s n ) g V3 n e s c J 
1 

< -. 

- 2 

Of course it is absurd. Now (13.91) yields 

Ho(X) = / ki(z, dy)k 2 (z, dx)iio(dz) 
Jn 



< / ki(z, dy)k 2 (z, dx)u,o(dz) + / ki(z,dy)k 2 (z,dx)u,o(dz) 

'Ai JA 2 

ki(z, dy)k 2 (z, dx)fj jQ (dz) 

A ± nA 2 

= fio(X)- [ (l-k 1 (z,{z}))(l-k 2 (z,{z})) t i (dz). 
Jx 

This equality asserts that there is Cl G 3B{X) with ^o(Cl c ) = so that ki(x,{x}) = 1 
or k 2 (x,{x}) = 1 holds for all x G Cl. Set Cl\ := {x G X ; ki(x, {x}) = 1}. Let 
l : I^IxIbe given by t(x) = (x,x). For A G ^(X), dSTTJ) yields 

no(A) = fi 2 (A n Xi) + /ii (A n x*) 



> / k 2 (z, {z})u, 2 (dz) + / ki(z,{z})fii(dz) 

^WAniiD+iiWinif)) 
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This estimate implies /i(t(f2 c )) = 0. Thus we have 

n(t(A)) = fi(t(A n Oi)) + n(t(A nn c x n n)) 



ki{z,{z})iJii{dz)+ I k 2 (z,{z})fx 2 (dz) 
AnQj JAnQfnn 

= ni{A n Oi) + /x 2 (A n Cii n n) 

>/i (^)- 

Thus we obtain fj,\ D = /x ° t -1 . □ 

Proof of Proposition 13.51 We use the same notation as in the proof of Lemma 13.61 
We denote p, = // — //o ° Note that p, is positive and it is absolutely continuous with 
respect to // by Lemma I3T61 In order to derive \i = fi f M , we consider the integration of ip s 
by /i' M for s G {s n } neN : 

ip s (x,y)u t M (dxdy) = / (p s (x, Rx)fi 1 (dx) - / ip s (x, Rx)fi (dx) 

\ x A ./ .V ./ .V 

<£> s (x, Rx)fj,(dxdy) — / <£>s(^, R%)fM) ° L^(dxdy) 

XxX JXxX 



<p s (x, Rx) n(dxdy) . 



XxX 



By virtue of (Al), we also obtain 



(p s (x,y)fi M (dxdy) = / (p s (y, Ry)fx(dxdy). 

XxX JXxX 



By Lemma [3.41 



= / ip s (x,y)n t M (dxdy) - / tp s (x,y)u.(dxdy) 

JXxX JXxX 



= ol {<p s {x,Rx) + <p s (y,Ry) -2<p s (x,y)}p(dxdy). (3.10) 

^ J XxX 

By (Al) and (j3TT|l . x X) = //(X x X ± ) = 0. This fact together with flBTTUD yields 

{y? s (x, Rx) + ify) - 2ip s (x, y)}fl(dxdy) = 0. (3.11) 



Let be as given m Lemma EH For zeliH 1^, 

ife) = P z [Z(s) G X 1 ] - F Rz [Z(s) G XJ 
= ¥ z [Z(s) e Xt] - ¥ z [Z(s) e X 2 }. 

Thus, for x,y £ (X x x X x c ) n (X w x RX®), 

(p s (x, Rx) + (p„(y, Ry) - 2cp s (x, y) 

= F x [Z(s) G Xi] - P„[Z(s) G Xi] + P„[Z(s) G X 2 ] - F x [Z(s) G X 2 ] 
-2 sup |P x [Z(s) G A]-P„[£(s) G A] | 

< 0. (3.12) 
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Note that 

x RX^ f ) < //i((l (t) ) c ) + fx 2 ((RX®) c ) = 

since fi E <^(/ii,/i 2 ). Hence there is E s C {X x x Xf) n (XW x RX®) with /t(£ s c ) = so 
that the equality holds in f[3TT2D for (x,y) E E s . Let E = E® n (fLew^n)- Here sW is 
given in Definition 12.51 associated with P G ^(P X1 , P X2 ). Then jl(E c ) = and 

~ ||F X o Z( Sn )- 1 - Pj, o Z( S „)- 1 || var = ¥ x [Z(s n ) E X x ] - F y [Z(s n ) G Xi] 

= Pj,[Z(s n ) G X 2 ] - P x [Z(s n ) G X 2 ] 

for all (a;, y) E E and n G N. Hence the property of immediately implies x = Ry for 
every (x, y) E E (cf. Remark |2~771 (i)). It yields fi = fi f M . □ 

Proof of Theorem [Ml Let t E (0, oo] and P G ^(P^P^) n %(¥ Xl1 V X2 ) maximal at 
each t E (0, to)- Note that 

T{Z U Z 2 ) = inf{s > ; Z 1 {s) = Z 2 (s)} P-a.s. (3.13) 

holds since P is maximal. Take s, t > with s + 1 < to- By the maximality of P at s + t, 

<p a +t(xi,x 2 ) =F[T{Z U Z 2 ) >s + t] 

= E[F[T(Z 1 ,Z 2 ) >s + t\^*}}. 

Since P G %{F X1 ,F X2 ), (Q yields F[T(Z l ,Z 2 ) > s + t\3? t *} > <p s (Z 1 (t),Z 2 (t)). In 
addition, by (13.31) . E[ip s (Zi(t), Z 2 (t))] > <p s+t (xi,x 2 ). Hence we obtain 

E[ip s (Z 1 (t),Z 2 (t))] = if s+t (x h x 2 ). 

Letting s — > 0, Proposition 13.51 yields Po (Zi(t), Z 2 (t)) _1 = /i^. Since £ G (0,t ) is 
arbitrary, it implies that 

F[Z 2 {t) = Z 1 (t) or Z 2 (t) = RZ 1 (t) for all t E (0,t )] = 1. 

Recall that r is the first hitting time of Z\ to if. The above equality implies that r equals 
the first hitting time of Z 2 to H P-almost surely. In addition, by (A2), for each t E (0, to), 

{t < t} c {Z 2 (t) = RZx{t)} P-a.s.. (3.14) 

Thus it suffices to show that 

{t < t} C {Z 2 (t) = Z 1 (t)} P-a.s.. (3.15) 

Note that implies P[Z x (t) = Z 2 (t)] = a4(X). By the maximality of P, Lemma O 
and (13TH . 

P[T(Z 1; Z 2 )<t] = l- ip t ( Xl ,x 2 ) = 1 - F Xl [Z(t) E X x ] + P, 2 [Z(t) G X x ] 

= P, 1 [Z(f)eX 1 c ] + P I3 [Z(t)el 1 ] 
= ^(X). 
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It means 



F[Z 1 {t) = Z 2 (t)]=F[T{Z 1 ,Z 2 )<t]. 



(3.16) 



Since we have ( 13.141) 



{Zi(t) 



-- Z 2 (t)} c {t < t} P-a.s., 
{r < t} C {T(Zi, Z 2 ) < i} P-a.s.. 



The second inclusion follows from (13.131) . Combining them with (I3.16p . we obtain (13. 1 5)) 



4 Examples and counterexamples 

Let us consider several examples of A and Z with (Al) and (A2) for given x%,x 2 G A. 
First we state a sufficient condition for (12.91) to be satisfied. A key ingredient is the 
Varadhan type short time asymptotic behavior of transition probabilities (14.11) . In order 
to state it in a general form, we introduce some terms concerning the metric geometry. 
Let (X, d) be a metric space. We call a curve 7 : [0, 1] — > X geodesic if, for each t G [0, 1], 
there exists 5 > so that d(y(t), 7(s)) = \t — s|d(7(0), 7(1)) for \t — s\ < 5. Recall that 
a metric space (A, d) is geodesic when, for each x, y G A, there is a rectifiable curve in 
A whose length realizes d(x,y). Note that such a curve always becomes a geodesic by 
a suitable re-parameterization. We call it a minimal geodesic joining 7(0) and 7(1). A 
geodesic metric space (A, d) is said to be non-branching when, for any two geodesies 7 
and 7' of the same length with 7(0) = 7'(0), we have inf {t > 0; j(t) 7^ j'(t)} = or 00. 
Here we follow the usual manner inf = 00. 

Theorem 4.1 Let (A, d) be a non-branching geodesic metric space and (Z, P x ) a diffusion 
process on it. Suppose that (Al) and (A2) hold for given Xi,x 2 G A. In addition, we 
assume the following for each t > 0: 

(i) the support of the law of Z(t) under F X1 equals X, 

(ii) there exist 

• an increasing function p : (0, 00) — > (0, 00) with lim s _>o p(s) = ; 

• a strictly increasing function \I> : [0, 00) — > [0, 00), 



and it completes the proof. 



□ 



• a sequence {s n } ne ^ of positive numbers with lim^oo s n = 

• G 38{X) with BY® = Y® and ¥ Xl [Z(t) (£ Y®] = 



and 



so that 



lim p(s n ) logP x [Z(s n ) G A] = %(d(x, A)) for each A G 3§{X) 



(4.1) 



n 



■00 



holds for x G Y® . 



Then holds. 
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Proof. Take z G X^nY^ for t > 0, where X® is as in LemmaEIB Then P 2 [Z(s) G A] = 
Wrz [Z{s) G RA] holds for each A G 3&{X). Thus flU} yields A) = d(Rz, RA). By 
taking A = B r (w), a ball of radius r > centered at w G X, and taking r — > 0, we obtain 
d(z,w) = d(Rz, Rw). Since R is continuous, the condition (i) implies that R acts on X 
as isometry. Let x G n X 1 and y G n X 2 . Suppose ([22D and (|S7SD holds. Take 
z G X\. Since Xi is open, B r (z) C Xi for all sufficiently small r > 0. For such r, ( 14. ip 
and ( 12. 7p implies that, 

*(d{x,B r {z))) = - lim p{s n )log¥ x [Z{s n ) G B r (*)] 
< - lim p{s n )logF y [Z{s n ) G B r (*)] 

n— >oo 

= *(d(j/,fl P (;z))). (4.2) 

It immediately implies 5 r (z)) < d(y, B r (z)). By taking r — > 0, we obtain d(x,z) < 
d(y,z). In the same way, for z G X 2 , we obtain c?(x, 2) > d(y,z). These two estimates 
yield 

d(x, z) = d(y, z) for z G H. (4.3) 

Take a minimal geodesic 70 : [0, 1] — > X joining x and y. By the remark after (A2), 
there exists to G [0,1] so that 70 (to) £ H. By ( 14.31) . we have d(x,jo(to)) = d(y, 7o(to))- 
In addition, t — 1/2 follows. Again by (14.31) . d(x, j (l/2)) = d(x,H) holds. Let 71 be a 
curve joining x and Rx given by 



7i (t) 



7o (t) if t e [0,1/2], 

i2( To (l-t)) if te (1/2,1]. 



Then 7x is a minimal geodesic joining x and -Rx because d(x, 7i(l/2)) = d(Rx, 7i(l/2)) = 
d(x,H). Since X is non-branching, we obtain Rx = y. Thus, once we set = (y® fl 
X x ) x (yW n X 2 ), P[Z*(t) G (SW) c n Xi x X 2 ] = holds for each P G ^(P^, P l2 ). This 
means tf{F Xl ,F X2 ) = tf(F xl , F X2 ) and therefore O holds. □ 

Remark 4.2 If our diffusion process Z(t) has a continuous transition density p t (x,y) 
with respect to a Radon measure m, that is, ¥ x [Z(t) E A] = j A p t (x,y)m(dy), Then ( 14. ip 
in Theorem 14.11 is replaced as follows: 

- lim p(s n ) \ogp Sn (x, y) = $(d(x, y)) for each y G X. (4.4) 

n— >oo 

Indeed, (12.71) and (12TS1) imply that p Sn (x, z) > p Sn (y, 2;) for z G X : and p Sn (x, z) < p Sn (y, z) 
for z G X 2 . Thus the same proof works. 

Corollary 4.3 Let X be a complete Riemannian manifold and Z(t) the Brownian motion 
on X. Assume X to satisfy (Al) and (A2). Then (12. 9p follows. 

Proof. In this case, Z has a continuous transition density p t (x,y). Letting p(s) := 2s, 
*ff(v) = v 2 and any sequence {s n } ne jj with lim^oo s n = 0, (14.41) follows from [T5] for every 
x G X (see Remark 14.21) . The condition (i) in Theorem 14.11 comes from strict positivity of 
the transition density. It is well-known that all properties imposed on X in Theorem 14.11 
hold. □ 
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We can also apply Theorem 14.11 to Alexandrov spaces. These metric spaces are an gener- 
alization of a complete Riemannian manifold with sectional curvature bounded below (see 
for details). 

Corollary 4.4 Let (A, d) be an Alexandrov space and Z(t) a diffusion process on X 
corresponding to a canonical regular Dirichlet form on X constructed in [7^/ (see [73] / 
also). Assume X to satisfy (Al) and (A2). Then (12. 9p follows. 

Proof. In this case, there is a continuous transition density p t (x,y) of Z. Letting 
p(s) = 2s, = v 2 and any sequence {s n } ne N with lim^oo s n = 0, (14. 4p follows 

from Corollary 2 of [16J for every i6l (see Remark 14. 2p . As in the case of Riemannian 
manifolds, the condition (i) follows from positivity of the transition density. By definition, 
X is a geodesic space. The curvature condition on X easily implies that X is non- 
branching. Thus we can apply Theorem 14.11 □ 

Remark 4.5 Let A be a Riemannian manifold and Z the Brownian motion on it. (i) 
If (Al) and (A2) are satisfied for given initial points, then the argument in the proof 
of Theorem 14.11 implies that R is isometry. In this case, if is a totally geodesic smooth 
submanifold of A (see [10J p. 61, for example). In particular, H becomes a complete 
Riemannian manifold. In addition, H is of codimension 1. (ii) If (A2) are satisfied with 
respect to an isometry R with R o R = id, then (Al) follows for each X\,X2 G A with 

RX\ = X2- 

In what follows, we will see some manifolds satisfying the conditions (Al) and (A2). In 
all cases, we assume Z to be the Brownian motion. 

Example 4.6 We consider the case A is an irreducible Riemannian global symmetric 
space of constant curvature. We will review that (Al) and (A2) are satisfied for every 
distinct pair x\, X2 G A of starting points in these cases. By Remark 14.51 (ii). It suffices to 
find an isometry R with RoR = id, Rx% = X2 satisfying (A2). The flat case, i.e. A = M n , 
is considered in [6]. 

In the case of positive curvature, A is a sphere: 

X = S n = {z = (z , Zl ,..., z n ) G R n+1 ; z 2 + ■ ■ ■ + z\ = r] 

with a metric induced from the canonical metric on M n+1 . Take xi,x% G A with x\ ^ x%- 
Then we can easily verify that the restriction of the reflection in M n+1 with respect to a 
hyperplane fulfills all of our requirements. 

In the case of negative curvature, A is a hyperbolic space: 

X = U n = {z = (zq, zx,..., z n ) G R n+1 ; -z 2 + z\ + • • • + z 2 n = -r, z > 0} 

with a metric induced from the Lorentz metric on M n+1 . Take X\,X2 G A with X\ ^ X2- 
Let m be the midpoint of x\ and X2- By homogeneity, we may assume m — (r, 0, . . . , 0). 
By arranging the chart appropriately, we may assume x\ = (zq, z\, 0, . . . , 0). Then X2 = 
(zo, —zi, 0, . . . , 0). Set R : (zq, z\, . . . , z n ) t— > (z , —z\, z 2 , . . . , z n ). Then R fulfills all of our 
requirements. 
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Remark 4.7 The converse of Example 14.61 is true for an irreducible Riemannian global 
symmetric space A in the following sense. If there exists an isometry R satisfying (Al) 
and (A2) for some pair Xi,x 2 G A, then X must be of a constant curvature. It follows 
from the result in [7] (cf. Remark |4.5f i)). 



Example 4.8 Let us consider 2-dimensional torus X = T 2 = (M/~) 2 , where ~ identifies 
x with x + n for each x G K and jiGN. Let 7r : M 2 — > T 2 be the canonical projection. 
We denote 7r(a?) by [s]. Take Xi,x 2 G A with xi 7^ X2. By arranging an appropriate 
chart, we may assume that x\ = [(a, 0)] and x 2 = [(0,6)] for < 6 < a < 1/2. Let 
K = {z G T 2 ; Xi) = d(z,x 2 )}. If (Al) and (A2) are satisfied for Xi,x 2 G A, then i? 
is isometry and H C K must hold (cf. Remark 14.51) . In the following, we will write K 
explicitly. First we consider the case 6 7^ 0. Take six points z\, z 2 , z$, Z4, z$, z% G M 2 as 
follows: 




Let Uj be a line segment in M 2 whose endpoints are Z{ and Zj. Then K = vr(/i2U/23U/45U/56) 
holds. We can easily verify that K has singularity at [z 2 ] or [z 5 ] (see Fig.l) and H cannot be 
contained in K by Remark l4.5( i). Thus there is no reflection structure. Next we consider 
the case 6 = 0. Then we have K = vr({(a/2,g) ; q G [0,1]} U {((1 +a)/2,q) ; q G [0,1]}) 
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(see Fig.l). Thus a map R : T 2 — ► T 2 denned by g)]) = [(a — p, g)] satisfies (Al) 

and (A2) with H = K. 

Example 4.9 Let A be a complete Riemannian manifold given by the direct product 
of two manifolds Y\ and Y 2 . We assume that the Riemannian metric on X has a form 
h(y 2 )cri(dyi ® dyi) + a 2 (dy 2 <8> cfa/2), where /i is a positive function on Y 2 and Oi is a 
Riemannian metric on Y{. We also assume (Al) and (A2) on (Y\,ai) for given starting 
points jji, y^ G Vi. Then we can extend the reflection structure on Y\ to A in a natural 
way. As the result, for any y^ G Y 2 , (Al) and (A2) are satisfied for (y^\ y^), (y^\y^) £ 



The same argument works for rotationally symmetric manifolds. Take a function 
h : (0, 00) — > (0, 00) which has a smooth extension to [0, 00) and satisfies h(0) = and 
h'(0) = 1. By using h, we define a metric on (0, 00) x S n_1 by ds 2 = dr 2 + h{r)d9 2 for 
(r, 6) G (0, 00) x § n_1 . Let X be the completion of (0, 00) x S n_1 by adding one point o. 
o is the limit of (r, 6) as r — ► for each 9 G S n_1 . Take x\ = (r, 81) and x 2 = (r, 9 2 ) for 
some r G (0, 00) and 9%, 9 2 G S n_1 . Then, the Brownian motion with starting points X\ 
and x 2 satisfies (Al) and (A2). 

Next we give two simple examples satisfying (Al) and (A2) while the uniqueness of 
maximal Markovian coupling does not hold. 

Example 4.10 (Fig.2) Let Y x and Y 2 be two copies of T 1 = [0, l]/~. For ^ = 1/26^ 
(i = 1,2), we set X = Yi U Y 2 j ~ where ~ means the identification of z\ and z 2 . Set 
x\ — G Y\ and x 2 = G Y 2 . We identify a function / on X with a function / on Y\ U Y2 
with /(^i) = f(z 2 ). We define a bilinear form <^ on < / ; / is smooth on Y\ U 1^ \ by 



We can easily verify that $ is closable in L 2 (X, dx) and its closure defines a Dirichlet form 
on X. Thus we can define the corresponding diffusion process ({Z(t)} t >o, {P x } x( z X ) ( see 
[1]). By using identity maps i\ : Yi — > F2 and 62 : 5^ — > li, we define : A — > A by 



A. 
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Rx = Li(x) if x e Yi. By definition of R and F Xi (i = 1,2), (Al) and (A2) holds. In this 
case, H = {z{\ = {z 2 }. Let us define a map n : Y 2 — > F 2 by ^(x) = 1 — x. We define a 
Markovian coupling P £ ^(P I1 , P X2 ) as the law of (Zi, Z 2 ) where Z\ is a copy of (Z, F X1 ) 
and 

I -Zl(t) if t > T. 

Then clearly P ^ F M and 

P [T(Z U Z 2 ) >t]= F X1 [r>t]= P M [T(Z U Z 2 ) > t] 

for each t > 0. Thus P is also a maximal Markovian coupling. Note that rj o R also 
satisfies (Al) and (A2) instead of R in this case. 

Example 4.11 (Fig. 3) Next example is a tree. The space X, given in Fig. 3, is a union 
of nine copies of the unit interval [0, 1] with some identification of these endpoints. X is 
naturally regarded as a metric space. As in Example I4.10[ we can construct a canonical 
Dirichlet form and the corresponding diffusion process on X. Let X\ = pu and x 2 = p 22 . 
There is an isometry R : X — > X so that R(pn) = p 22 , R(pn) = P21, R(Pi) = P2 
and R fixes all other endpoints. Then (Al) and (A2) holds. Let n be an isometry 
so that f](p 2 i) = p 22 and 77 fixes all other endpoints. We define a Markovian coupling 
P G ^(PjpP^) as the law of (Zi, Z 2 ) where Z\ is a copy of (Z,¥ Xl ) and 

(R(Z 1 (t)) if t<r {pih 
Z 2 (t) =Lo R{Z 1 {t)) if r {pi} < t < r {po} , (4.6) 
(Zx(t) ift>r {po} , 

where ts x \ is the first hitting time to x. Then clearly P 7^ Fm and 

P [T(Z 1 , Z 2 ) >t]= F X1 [t {po} >t]= F M [T(Z 1; Z 2 ) > t] 

for each t > 0. Thus P is also a maximal Markovian coupling. Different from Exam- 
ple 14.101 this example essentially has only one reflection structure. 

These examples reveal that maximal Markovian coupling may not be unique if the 
underlying space is more singular than Riemannian manifolds or Alexandrov spaces. One 
characteristic property which is common to those examples is the existence of branching 
geodesies. But, in general, non-branching property of geodesies is not necessary for the 
uniqueness of maximal Markovian coupling. To see this fact, we consider the Brownian 
motion on 2-dimensional Sierpinski gasket. 

Take three points pi,P2,P3 £ P 2 with \pi — Pj\ = 1 for all i 7^ j. Let us define a 
contraction map : M 2 — > M 2 for i = 1, 2, 3 given by ^i(x) = (x — Pi)/2 + Pi. Obviously, 
Pi is the unique fixed point of The Sierpinski gasket is a unique compact set in R 2 
satisfying X = Ui=i ^iPO ( see Fig-4). For detailed properties of the Sierpinski gasket, 
see [9] for instance. We set Vq = {^1,^2,^3} and V n = UiLi ^i(^n-i)- The Brownian 
motion ({Z(t)} t >o, {F x } x( zx) on X is given by a suitable scaling limit of a continuous 
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a geodesic joining p 3 and 



* 2 (Pl) 



Fig.4 



time random walk on as n — > oo (see [U E]). There is a reflection i? on R 2 so that 
R{p\) = P2- We denote the fixed points of R by H. The map naturally induces a 
reflection on X so that its fixed points H coincides with X f] H. Moreover, X and 
({Z(t)}t>o, {IPijxex) fulfills (Al) and (A2) for X\ = p\ and x 2 = P2- As shown in [9], 
there is a unique distance donX, called shortest path metric, such that it satisfies 

(i) (X, d) becomes a geodesic metric space, 

(ii) d(pi,pj) = 1 for each i 7^ j, 

(iii) d(zi, z 2 ) = 2d(V i (z l ), ^i{z 2 )) for z x , z 2 e X and i = 1, 2, 3. 

Theorem 4.12 Let X be the Sierpinski gasket as defined above. Then ( 12. 91) holds. 

Proof. Let p t (x,y) be the transition density of the Brownian motion. Then, the main 
theorem of [IT] asserts that, for each u > and 1,1/6 X, 



lim 

n— >oo 



m) logp ( 2/5)»u(x,y) = d(x,y) dw/{d ™ l) F 



u 



d(x,y)J ' 



(4.7) 



where d w > 2 is the walk dimension of the Sierpinski gasket and F is an implicitly 
determined, non-constant, positive continuous function on (0, 00). For our aim, we need 
a refined observation on F. By the definition of F in [TT] . 

(4.8) 



F(v) = v 1/{dw - l) sup {K(s) - vs} 



s>0 



for some positive, concave and real analytic function K(s) on (0, 00). Thus, 



d{x,y) dwl ^- x) F 



u 



d(x,y) 



u 



l/(d w -l) 



sup {d(x, y)K(s) — us} . 

s>0 
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Since F is continuous on (0,oo), there is s v G (0, oo) for each v G (0, oo) so that 
K(s„) — vs v — sup s>0 {K(s) — vs} holds. Indeed, if there exists a sequence {s n }neN 



with linin^oo s n = oo so that 

lim (K(s n ) - vs n ) = sup {K(s) - vs} , 

S n ^OO S> Q 

then F(v') = oo for every v' < v. These observations imply that, for < a < b, 

sup {aK(s) - us} = aK(s u/a ) - us u/a 

< bK(s u/a ) - us u/a 

< sup {bK(s) — us} . 

s>0 

It means that the right hand side in (j4.7p is strictly increasing with respect to d(x,y). 
Thus, the same argument as given in Theorem 14.11 yields that x G X% and y G X 2 with 
(I2TD and (TjLBD satisfies 

d(x, z) = d(y, z) for all z G H. (4.9) 

To complete the proof, we show that (14.91) implies x = Ry. It suffices to show that w = w' 
holds when w, w' G X 2 with d(w, z) = d(w f , z) for z = p 3 or z = \E , 2 (Pi)- in this case, we 
have 

we^ 2 (X) & d(w,p 3 )>l/2,d(w,y 2 (pi))<l/2, 
we* 3 {X) & d(w,p 3 )<l/2,d(w,^ 2 (p 1 ))>l/2. 

Thus, w,w' G ty 2 (X) or w,w' G *ff 3 (X). In particular, w = w' = ^f 2 (p 3 ) if and only if 
d(w,p 3 ) = d(w, ty 2 (Pi))- Now we assume w G ty 2 (X)\ty 3 (X). To see the argument below, 
we easily find that the same argument also works for the case w,w' G ^ 3 (X) \ ^ 2 {X). 
Since (X, d) is a geodesic space, d(w,p 3 ) = d(w, ty 2 (p 3 ))+l/2 and therefore d(w, ^ 2 (p 3 )) = 
d(w f , ^ 2 (p 3 )). Then we have 

we^ 2 o^ 1 (X) d(tu, * 2 (pi)) < l/4,d(tw,* 2 (p3)) > V 4 > 
w e^ 2 oty 2 (X) & d(w,y 2 (pi)) >l/4,d(w,y 2 (p 3 )) >l/4, 
we^ 2 o^ 3 (X) d(w,^ 2 (pi)) > l/4,d(w,y 2 (p 3 )) < 1/4. 

Thus w, w' G ^ 2 (^i(X)) for some i G {1, 2, 3}. In particular, w — w' when w G V 2 . Since 
we have 

d(w, * 2 (pi)) = *2 o *i(p0) + rf(^ 2 O *i(pi), # 2 (Pl)) S 
* 2 (pa)) = d(tw, ^2 o *i(p 3 )) + rf(^ 2 o *<(p 3 ), *a(ps)) 

when w G \l/2 ^i(X), the same argument as above works by replacing ^ 2 (pi), ^ 2 (p 3 ) and 
y 2 (X) by ^ 2 o Vi(pi), ^2 o *i(p 3 ) and ^ 2 ° respectively. When w G UneN V n, such 

a recursive argument ends in a finite step with resulting w = w'. When w IJneN V n , we 
obtain a sequence {i n } n£ N with i n G {1, 2, 3} so that w, w' G ^ o \l> i2 o • • • o ^> in (X) for 
each n G N. Since f\eN ' ' ' ^i™P0 i s J us t one point, w = w' follows. □ 
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Remark 4.13 As shown in the above proof, (14.71) means that (14.41) holds with p(s) = 
s i/(d w -i)^ g ^ _ u ^2/5) n and ^{y ) = v dw ^ dw ~^ F(u/v). Thus the Sierpinski gasket satisfies 
all assumption in Theorem 14. II except for being non-branching. For example, we consider 
two minimal geodesies 71 and 7 2 . 71 joins p 3 and p 2 . 72 joins p 3 and ^(Pi) via ^2(^3)- 
Then both of 71 and 72 contains the minimal geodesic joining p 3 and ^2(^3) • Thus 
Theorem 14.121 is not a direct consequence of Theorem 14. 1L 



5 Kendall-Cranston coupling 

Let X be a d- dimensional complete Riemannian manifold and ({Z(t)} t > , {F x } x( zx) the 
Brownian motion on it. In this framework, we construct a Kendall- Cranston coupling 
following the argument due to von Renesse [17] . As we will see, his argument is based on 
approximation by coupled geodesic random walks. 

Let D(X) = {(x,x) G X x X ; x G X}. For each (x,y) e X x X \ D(X), we choose 
a minimal geodesic ^ xy : [0, 1] — > X of constant speed with 7x y (0) = x and 7^(1) = y. 
Let H xy be the hyperplane in T y X of codimension 1 which is perpendicular to 7ry(l) and 

G H xy . For each v G T X X, take a parallel translation by Levi-Civita connection along 
7 xy to T y X and reflect the resulting vector with respect to H xy . In this way, we obtain a 
new vector w G T y X. We define a map m xy : T X X — > T y X by m xy v = w. Clearly m xy is 
isometry. Take a measurable section ip : X — > G(X) to the orthonormal frame bundle 
6{X). Let us define maps $j : IxI-4 G{X) for i — 1, 2 satisfying 

$i(x,y) G X (X) x.yeXxX, 

$ 2 (x, y) G V (X) x, y G X x X, 

$ 2 (z, = m iC j / $ 1 (a;, y)u (x, y) G X x X \ D(X),u G R d , 

<&i(x,x) = $2(^5^) = <p{ x ) x £ X. 

We can choose j xy so that (x, y) 1— > 7^ is measurable as a map from X x X \ D(X) to 
C 1 ([0, 1] — >■ X) and j xy is symmetric, i.e. j X y(t) = Jyx{l — t). Also we can choose 
to be measurable for i = 1,2. Take a sequence of random variables {£n}neN uniformly 
distributed on <i-dimensional unit disk. Let us define a coupled geodesic random walk 
Z e (n) = (Zf(n), on X x X with step size e > and starting point (x, y) G X x X 

inductively by 

^ £ (0) = Or,y), 

Z £ (n + 1) = (exp Zf(n) ( £ Vd + 2$i(£ £ (n))£„+i) , exp Z|(n) (ev'dT^ 2 (Z £ (n))£ n+1 )) • 

Let T A (t) be the Poisson process with intensity A > independent of {£ n }neN- Then the 
sequence of processes {Z k 1/2 (r fc (t))} fcgN is tight in the Skorokhod path space D([0, 00) — > 
X x X) and Zf ' (T k (t)) weakly converges to the Brownian motion on X as k — > 00 for 

1 = 1,2. Let Z(t) = (Zi(t) , Z 2 (t)) be a (subsequential) limit of {Z k 1/2 (r k (t))} keN . Let a 
be the first hitting time of Z to D(X). We set Z(t) by 



Z(f) 



Z(t) if t < <t, 

{Z X {t),Z x {t)) ift><7. 
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We call Z(t) a Kendall- Cranston coupling. This is indeed a coupling of two Brownian 
motions starting at x and y respectively. Our choice of is a bit different from that in 
[TT] . where £ n is uniformly distributed on the unit sphere. But it does not matter since 
the same argument works. 

Theorem 5.1 Assume (Al) and (A2) for x±,X2 G X. Then a Kendall- Cranston coupling 
of¥ Xl and ¥ X2 is the mirror coupling defined by R. In particular, the Kendall-Cranston 
coupling is unique in the sense that it is independent of the choice of subsequences of 
approximating geodesic random walks. As the result, the Kendall- Cranston coupling is the 
unique maximal Markovian coupling of F Xl and F X2 . 

Proof. As shown in the proof of Theorem 14.11 R is an isometry on X. For x G X, set 
y = Rx. We claim 

dR{u) = m xy u for u G T X X. (5.1) 

In order to complete the proof, it suffices to show (15.11) . Indeed, since the equality 
R(exp z (w)) = exp Rz (dR(w)) holds for z G X and w G T Z X, (15.11) implies 

Z|(n) = RZ{{n) for n < T(Zf, Zf). (5.2) 

Note that the coupled geodesic random walks never meet under (15.21) . That is, 

P[T(Z 1 £ ,Z 2 £ ) = oo] = l (5.3) 

for each e > 0. This fact is shown as follows: by (15. 2p . Z\(T(Z\, Zf )) G H must hold if 
T(Zf, Zf ) < oo. Let v z>e be the law of exp z (e!;i). Then v Z£ (H) = for each z G X and 
e > since if is a submanifold of codimension 1 as mentioned in Remark 14.51 It easily 
implies (15.31) . Once we obtain (15.21) and (15. 3p . the central limit theorem for geodesic 
random walks yields that the full sequence of {Z k 1 2 \rt (t))} ken weakly converges to the 
image of the Brownian motion by the map z i— > (z, Rz) as k — ► oo. Thus a Kendall- 
Cranston coupling is unique and identical to the mirror coupling. 
Set 

7 — Ixy f° r simplicity. First we show 

fl( 7 (t))= 7 (l-t). (5.4) 

By the symmetric choice of 7^, we may assume I := inf zg # d(x, z) < inf ze ^ d(y, z) without 
loss of generality. By (A2), / < 00 holds. Take to G [0,1] so that rf(x,7(to)) = I- Let 
71 : [0, 2t ] — > X be a curve joining x and y given by 



7i(s) 



f 7 (s) sG[0,t ], 
(R(j(2t -s)) se(t ,2t 



Then the length of 71 equals 21 and the minimality of 7 implies 21 = d(x, y) and to = 1/2. 
Moreover, 7x(l/2) = 7(1/2) must hold and therefore 71 = 7. It proves (15 .4p . Note that 
the above discussion implies 7 (£) G Xi for t G [0, 1/2) and 7 (i) G X 2 for t G (1/2, 1]. 
Next we show (15. ip in the case u = 7(0). It easily follows from (15 .4p . that is, 

dR(u) = dR(7(0)) = -7(1) = m xy u. (5.5) 
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Finally we prove (15. ip for u _L 7(0). Let // s ^ T^^X — > T^mX be the parallel translation 
along j\[s t t\- It suffices to show that 

//l,l/2 o dR o //i A o(/i) = /i (5.6) 
for each /i G T 7 (i/ 2 )X with /i _L 7(1/2). Indeed, once we prove it, 

dR(u) = dR(// i/ 2 fi O //qa/ 2 {u)) = /i/2,l / ,l/2(w) = /o,l( M ) = 

Now we show (15.61) . Take e > so that the exponential map exp 7 n i 2 \ : T 7 (i/ 2 )X — > X 
is diffeomorphic on 2e-ball centered at G T 7 (i/ 2 )X. We may assume that \h\ = e. Let 
h! = exp-J^ (7(1/2 - e/d(7(0), 7 (l)))). Note that // 1)1/2 o di? o // 1/2)0 (/i / ) = -/i'. We 
consider a curve c : [0, 1] — > X given by c(t) = exp(cos7rt h' + sirnrth). Since c(0) G Xi 
and c(l) G X 2 , (A2) yields that c intersects -£/\ By the choice of h', R(c(t)) ^ c(t) if 
t 7^ 1/2. Hence c(l/2) G i?. It implies ( 15.61) and therefore completes the proof. □ 

Acknowledgment. The problem treated in section [5] was suggested by K.-Th. Sturm. The 
author would like to thank him for his advice. 
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